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GIANT PLANET OCCURRENCE IN THE STELLAR MASS-METALLICITY PLANE 
John Asher Johnson^-^, Kimberly M. Aller^'\ Andrew W. Howard'"', Justin R. Crepp^ 

Draft version May 25, 2010 

ABSTRACT 

Correlations between stellar properties and the occurrence rate of exoplanets can be used to inform 
the target selection of future planet search efforts and provide valuable clues about the planet formation 
process. We analyze a sample of 1266 stars drawn from the California Planet Survey targets to 
determine the empirical functional form describing the likelihood of a star harboring a giant planet as 
a function of its mass and metallicity. Our stellar sample ranges from M dwarfs with masses as low as 
0.2 Mq to intermediate-mass subgiants with masses as high as 1.9 Mq. In agreement with previous 
studies, our sample exhibits a planet-metallicity correlation at all stellar masses; the fraction of stars 
that harbor giant planets scales as / c>c lO^-^'^®/^]. We can rule out a flat metallicity relationship 
among our evolved stars (at 98% confidence), which argues that the high metallicities of stars with 
planets is not likely due to convective envelope "pollution." Our data also rule out a constant planet 
occurrence rate for [Fe/H] < 0, indicating that giant planets continue to become rarer at sub-Solar 
metallicities. We also find that planet occurrence increases with stellar mass (/ oc Af^,), characterized 
by a rise from 3% around M dwarfs (0.5 M©) to 14% around A stars (2 M©), at Solar metallicity. We 
argue that the correlation between stellar properties and giant planet occurrence is strong supporting 
evidence of the core accretion model of planet formation. 
Subject headings: Methods: Statistical — Stars: Planetary Systems — Stars: Statistics 



1. INTRODUCTION 

Mass and chemical composition are key quantities in 
the formation, evolution and fate of stars. A star of 
a given age is, to first order, characterized by these 
two physical parameters, and the influences of mass and 
metallici ty extend to the formation and evolution of 
planets (|Johnsonll2009f ). Even the first handful of ex- 
oplanet discoveries revealed that the likelihood of a star 
harboring a planet was clos ely tied to ste llar iron con- 
tent, or metallicity [Fe/H] (Gon zaleall997 1. Subsequent 
studies of larger samples of stars using uniform spectro- 
scopic modeling techniques found that giant planet oc- 
currence increases sharply for stellar metallicity in excess 
of the Solar value, rising from 3% for [Fe/H] < to 25% 
for [F e/H] > +QA (jSantos et al.ll200llFischer fc Valentil 
[20051 hereafter FV05). 

In addition to informing models of p lanet formation 
(llda fc LinI 120051 : iMordasini et al.l 120091 : iJohansen et al.l 
l2009f ). theplanet-metallicitv correlation (PMC) has pro- 
vided a guide for the target selection of subsequent planet 
searches. The Next 2000 Stars (N2K) and Metallicity- 
Biased CORALIE surveys leveraged the higher metallic- 
ities of their samples to detect large numbers of close- 
in planets, many of which transit their host stars and 
thereby yield key ins ights into the interior structures 
of Jovian exoplanets (iFischer et afl 12 005: B ouchv et al.l 
120051: 1 Johnson et al.l [20061: iMoutouet al.l 120061) . Indeed. 



studies of known transiting planets have revealed ev- 
idence of a correlation between plan etary core mass 
and the met allicity of their host s tars (iSato et al.l 120051: 
Torres et al.] 120081: iGuillot et all 120061: [Burrows et al.l 
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While the first planet detections yielded a definitive 
correlation between giant planet occurrence and stellar 
metallicity, until recently very little was kn own about 
the effects of stellar mass ([Laws et al.l 120031 ) . The first 
Doppler-based planet surveys concentrated primarily on 
stars with masses similar to the Sun, both because 
it was desirable to find Solar System analogs and be- 
cause Sun-like stars make excellent planet-search targets. 
Compared to more massive stars, dwarfs with masses 
within 1.0 ± 0.2 M© are relatively numerous, have cool 
atmospheres and slow rotational velocities (Kot sin i < 
5 km s~^). The latter two features result in a high den- 
sity of narrow absorption lines in the spectra of Sun-like 
stars, which is ideal measuring for stellar Doppler shifts 
to high precision. 

Stars at the lower end of the mass scale (the K and 
M stars) are even more numerous than the Sun and 
they also display large number of narrow absorption fea- 
tures in their spectra. However, most low-mass stars are 
optically faint (Vmag ^ 10) and are thus not included 
in large numbers in most Doppler surveys. The faint- 
ness of late-K and M-t ype dwarfs can be overcome by 
using larger telescopes (jButler et al.ll200^ iBonfils et al.l 
'2005H) , and more re cently by observing at infrared wave- 
lengths (jBean et al.li2009 ') . ]3espite the small numbers of 
M dwarfs thus far monitored by Doppler surveys, one re- 
sult has become apparent: M dwarfs harbor Jovian plan- 
ets very infrequently. Only eight systems containing one 
or more giant planets have been found among the « 300 
M dwarfs on various Doppler p rograms (Johnso n et al.l 
l2010hl : [Hag"highiDour et aL[[20T(il) . 



It was originally thought that the paucity of Jupiter- 
mass planets around M dwarfs wa s due to a metallicit y 
bias among nearby, low- mas s stars (iBonfils ct al."2005aV 
However, a recent study bv iJohnson fc A dds (2009) re- 
vealed that M dwarfs likely have the same metallicity 
distribution as Sun-like stars, and stars with masses 
Mi, < 0.5 M0 are 2-4 times less likely than Sun- like stars 
to have a Jupiter ([Johnson et al.ll2007air2010b[ ). 

At the other end of the mass scale, the problems 
inherent to massive, early- type stars can be over- 
come by o bserving targ ets at a later stage o f their 
evolution (IHatzes et all 120031: ISetiawan et all 120051: 



Sato et all 120051: iReffert et all 120061: IJohnson et all 
2007bl: iN icdziclskTeran [2007t 



Nicdzicls ki et al.l 
Dollinger et al.i 120091) 



Once 



iLiuetalJ 120081: 

stars exhaust their 



core hydrogen fuel sources they move off of the main 
sequence, become cooler, and shed a la rge fraction of 
their primordial angular mom entum (jGrav fc Nagail 
119851: Ido Nascimento et al.|[2000t) . The effects of stellar 
evolution transform a 2 M© star from an A-type dwarf 
with V^otsini ~ 100 km 8"^ and Teg = 8200 K, to a 
K-type subgiant or giant w ith Vrot sin i < 2 km s~^ and 
Jeff ~ 4800 K jde Mede ir^os_eLalJ [19971: IGirardi et aLJ 
120021 : [Sandage et al1l2003l ). Surveys of "retired" massive 
stars have resulted in the discovery of w 30 Jupiter-mass 
pl anets with well-cha racterized orbits (see e.g. Table 1 
of lBowler et~aI1[20lol) . 

Usin g a sample of stars s panning a wide range of 
masses. TJohnson et al.l ([2007a) measured a positive corre- 
lation between stellar mass and the fractio n of stars with 
detectable planetfl In a related study, iBowler et al.l 
([20101) measured a planet occurrence rate of 26~t^% 
among a uniform sample of 31 massive subgiants. Fur- 
thermore, based in part on a study of planets around 
K giants in nearby open clusters, iLovis fc Mavor (2007) 
found that the average planet mass increases as a func- 
tion of stellar mass, indicating that gas giant planets 
become either more massive on average, or more nu- 
merous (or both) with increasing stellar mass (see also 
IBowler et al. 2010). 

The observed correlation between stellar mass and the 
occurrence of detectable planets, like the PMC before 
it, has added an important new variable to models of 
planet formation. While the Sun and Solar-mass stars 
serve as important benchmarks for understanding the 
formation of our own planetary system, successful, gen- 
eralized planet formation theories must now account for 
the effects of stellar mass, and presumably by exten- 



sion, disk mass (iLaughlin et al.l 



iKennedv fc Kenvonll2008t iCurriel 



200llIdafcLml[2005[: 



2009). 



While previous studies have uncovered the existence 
of a positive correlation between stellar mass and planet 
occurrence, it is important to understand the underly- 
ing functional form of the relationship. For example, it 
would be advantageous to know whether the correlation 
is purely linear, or if it can instead be better described as 
some other functional form. Besides informing theories of 
planet formation, an improved understanding of the rela- 
tionship between planet occurrence and stellar mass will 
also help guide the target selection of future surveys, and 
aid in the interpretation of results of current and future 

^ In that study, "detectable planets" were defined as having 
Mp sini > 0.8 Afjup and a < 2.5 AU. 



planet-search efforts. Just as some previous Doppler sur- 
veys biased their target selection toward high-metallicity 
stars to increase their yield, future direct-imaging, as- 
trometric and Doppler surveys may benefit from concen- 
trating on more massive stars. This strategy has paid 
off for one high-contrast imaging survey, resulting in the 
detection of three giant plan ets around the A5 dwarf 
HR8799 ([Marois et all 120081) . Another example of an 
imaged planet is Fomalhaut b, which is a giant planet 
(< 3.3 Mjup) orbiting just inside of a debris disk of an 
ASVstar "(Kalas et al. 2008; Chiang ct al. 2009). Even in 
the cases when surveys do not yield detections, proper in- 
terpretation of null results requir es knowledge of t he ex- 
pected number of detections (e.g. lNielsen fc Closi [2009). 

Ascertaining the underlying form of the dependence 
of giant planet occurrence on stellar mass requires a 
larger sample than used in p r evious studies. Since the 
publication of IJohnson et al. l ([2007al) a sample of 240 
new intermediate-mass subgiants have been added to 
the California Planet Surv ey (CPS) at Keck Observa- 
tory [johnson^eF^I ([2010aD . At the low-mass end, two 
new giant planets have been discovered among the CPS 
Keck sample of M dwarfs (Johnson ct al. 2010b). Im- 
provements in our ability to estimate the metallicities of 
M dwarfs and massive evolved stars have provided vi- 
tal information about how to properly isolate the effects 
of stellar mass from the known effects of stellar metal- 
licity. With these tools at hand, we are now poised to 
make an updated evaluation of the relationship between 
stellar mass and planet occurrence. 

Our paper is organized as follows. In §[2]we present the 
characteristics of our three primary samples, including 
low-mass M dwarfs from Keck observatory; "Sun-like" 
late-F, G and K (FGK) dwarfs from the main CPS sam- 
ple; and massive, evolved stars from the Lick and Keck 
subgiant surveys. In §[3|we examine the separate effects 
of mass and metallicity on planet occurrence. In §[3|we 
present our Bayesian inference technique of measuring 
correlations between planet occurrence and stellar char- 
acteristics and we provide the best-fitting parameters for 
the measured relationship in § [5l We compare our results 
with previous work in § [6l Finally, we summarize our key 
results and discuss our findings in the context of the cur- 
rent theoretical understanding of planet formation in §[7| 

2. SELECTION OF STARS AND PLANETS 

Our goal is to measure planet occurrence as a func- 
tion of stellar properties. Care must be exercised in se- 
lecting the sample of target stars such that planets of 
a given mass and orbital semimajor axis could be uni- 
formly detected over the entire sample. The criteria for 
planet mass and semimajor axis translate into limits on 
velocity amplitudes, K, and orbital periods that can be 
tallied among the sample of planet detections. These 
criteria must be selected to ensure reasonably uniform 
detection characteristics across several Doppler surveys, 
which have different detection sensitivities and time base- 
lines. 

In what follows, we describe our selection of stars and 
planets from among the various CPS planet search pro- 
grams. The CPS is a collection of Doppler surveys car- 
ried out primarily at the Lick and Keck Observatories. 
The CPS target lists provide a large stellar sample with 
a wide range of masses and metallicities. Specifically, 
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Fig. 1. — H-R diagram showing our stellar sample (filled circles) 
and planet host stars (open circles). The solid line is the poly- 
nomial relationship describing the mean Hipparcos main sequence 
IWright) I I2004D . The excess scatter seen about the lower main se- 
quence B — V > 1.2 is primarily due to errors in the published 
V-band magnitudes and/or parallaxes for those faint stars. 

our stars lie in the ranges 0.2 < Mi,/MQ < 2.0 and 
— 1.0 < [Fe/H] < -1-0.55. The long time baselines ranging 
from 3 to 10 years, and Doppler precision ranging from 
1-5 m s~^ have resulted in a diverse and fairly com- 
plete sample of giant planets that have b een compiled 
in the C atalog of Nearby Exoplanets (^CNE. lButler et aLl 
(|2006bD l as updated by Wright et al. 2010 (in prep) in 
the Exoplanet Orbit Databas^^. 

2.1. Stellar Sample 

The low-mass stars in our sample are drawn from 
from the CPS Keck surve y of late-K and M-type dwarfs 
(|Rauscher fc Marcvl [20061 : IJohnson eFall I2010bn . This 
sample comprises stars with M^ < 0. 6 Mf^^ as esti- 
mated with the photometric calibration of lDelfosse et al.l 
()2000[) . We estimate the metallici ties with the broad- 
band photometric calibration of Johnson fc Appj ()2009D . 
which relates the metallicity of a star to its "height" 
(AMk) above the mean main-sequence in the {V — Ks, 
Mks} plane. 

The bulk of our Solar-mass F, G and K dwarfs are 
taken from the Spectroscopic Propertie s of Cool Stars 
catalog (SPOCS: rValenti fc Fischeil[2005l) . Most of these 
stars have masses in the range 0.8 < Mi^/M^ < 1.2. 
However, the SPOCS catalog contains some higher mass 
subgiants, which we fold into our high-mass stellar sam- 
ple described herein. The spectroscopic properties listed 
in the SPOCS catalog were measured using the LTE 
spectral synthesis software package Spectroscopy Made 
Easy ('SME: IVa lenti fc Piskunov 1996), as described by 
iValenti fc Fischer, f,2005i' ) and FV05. S tellar masses 
for th e SPOCS catalog are cataloged by iTakeda et al.l 
(|2007D , who associate the spectroscopic stellar properties 
to isochrones computed using the Yale Stellar Evolution 
Code ( YREC [AneFani2007l) . 

We select our high-mass stellar sample from the Lick 
and Keck Subgiant Planet Surveys. The sample selection 

^ [iittp : //exoplanets . org/ 1 



is descr ibed in IJohnson et al.l ()2006l ) and IJohnson et al.l 
(|2010al) . The masses and metallicities of the subgiants 
in our sample are estimated using SME and are listed in 
the fourth contribution to the SPOCS catalog (Johnson 
et al. 2010c, submitted). The majority of our subgiants 
have masses in the range 1.3-2.0 M©, with a tail in the 
distribution extending to 1.0 Mq. The metallicities of 
the subgiants range from [Fe/H] — —0.2 to -1-0.5. 

Our full stellar sample contains 1194 stars: 142 M 
and late-K dwa rfs from the Keck M Dwarf Survey 
(iButler et al.l l2006ai. 807 dwarf and subgiant stars from 



the original SPOCS catalog, and 246 subgiants from the 
SPOCS IV. catalog. Figure [T] shows our stars in the 
{Afy, B — V} H-R diagram. The open symbols are the 
positions of all of our stars, and the filled symbols are 
the stars known to harbor at least one detectable (giant) 
planet, as described in the following section. 

2.2. Planet Detections 

Following FV05, we restrict our analysis to systems 
with at least one "uniformly detectable planet," which 
we define as those with velocity semiamplitudes K > 
20 m s~^ and semimajor axes a < 2.5 AU. We de- 
creased threshold in K from the value used by FV05 
{K > 30 m. s~^) because of the increased Doppler pre- 
cisi on of HIRES since th e 2005 detector upgrade (see 
e.g. iHoward et al.l [20T0al) . For reference, at 1 AU and 
for circular orbits, semiamplitudes K = 20 m s~^ cor- 
responds to minimum planet masses Mpsini/ Afjup = 
{0.44, 0.82, 1.12} for M^/M© = {0.4, 1.0, 1.6}. 

Due to the limited time baselines of the Doppler sur- 
veys from which our targets are drawn, we also restrict 
our analysis to planets with a < 2.5 AU. This cri- 
terion is set primarily by our sample of intermediate- 
mass subgiants, which are on surveys with time base- 
lines ranging from 3-6 years. These criteria will, for 
most stellar masses, represent conservative cuts on the 
total number of giant planet detections. We defer the 
analysis of the frequency of less massive planets with 
Afpsini < 1.0 M.T„n or orbit s wider t han 2.5 AU to 
other studies (e.g. ISipusa et al][2008: Howard et al.ll2dol 



IGould et al.ll20lliCumming et al.ll2008D . 

To further ensure uniform detectability within our stel- 
lar sample we restrict our analysis to stars with a min- 
imum number of observations. For the low-mass and 
Solar-mass samples we require > 10 observations. For 
the high-mass subgiants we require > 6 observations; a 
smaller number owing primarily to the shorter time base- 
line of our Keck survey. We also require minimum obser- 
vational time baselines corresponding to our semimajor 
axis limit of a < 2.5 AU. Thus, for the M dwarfs we 
require a baseline > 6.3 years, using an average stellar 
mass M* = 0.4 Mq; > 4 years for the Solar-mass stars; 
and > 3 years for the subgiants with an average stellar 
mass of 1.6 M©. 

We compiled our sample of planet detections by cross- 
correlating our stellar samples with the Exoplanets Data 
Explorer, and recent planet announcements fr om the 
CPS (How ard et al.ll2010ht [lohnson et al.ll20T0allbl ). We 
augmented this list of secure detections with unpub- 
lished detections from the Keck Subgiants Planet Survey. 
These unpublished candidates all have more than 10 ob- 
servations over w 3 years, but lack strong enough con- 
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Fig. 2. — A plot of stellar mass (Af*) and metallicity ([Fe/H]) for the full stellar sample, comprised of 1194 stars (black dots), 115 
of which harbor at least one detectable planet (red diamonds). For visualization, we have divided the stellar sample into three broad 
groups: M dwarfs, FGK dwarfs, and massive "retired" A stars. The fraction / of stars with planets is printed above each group. The 
thick (black) lines each stellar mass group represent the best-fitting linear relationships between mass and metallicity (for the M dwarfs 
the lines represent the metallicity). The dashed (red) line is the best-fitting linear relationship between mass and metallicity for the stars 
with planets. For the M dwarfs we simply report the average metallicity for each population. The (blue) 2-dimensional error bars represent 
the typical measurement uncertainties. In each mass group, there is a systematic metallicity offset between the stars with and without 
planets. Discontinuities between the samples are not entirely physical, and are in large part due to the different target-selection criteria for 
the three surveys. 



straints on the orbital parameters for publication. How- 
ever, since the present study is concerned with planet 
occurrence we feel confident in including these secure, 
yet unpublished detections in our sample. All of the un- 
published candidates have radial velocity variations con- 
sistent with Doppler amplitudes and periods that meet 
our criteria for uniform detectability. 

Our sample of planet detections comprises 5 planets 
around M dwarfs, 74 planets around the SPOCS sample 
of FGK dwarfs, and 36 planets around subgiants. 

3. DISENTANGLING MASS AND METALLICITY 

In our analysis we treat stellar mass and metallicity as 
separate independent variables affecting the likelihood 
that a star harbors a planet. The validity of this premise 
rests in part on the analysis of FV05, who noted an ar- 
tificial correlation between mass and metallicity in the 
SPOCS sample that is due to the color and magnitude 
cuts used in the target selection: the more massive stars 
in the SPOCS sa mple have higher metallicities than the 
lower-mass stars (jSantos et al.l 120041 : iMarcv et al.l 120051 
;FV05). This selection effect is clearly seen in our up- 
dated data set shown in the middle panel of Figure [2l 



However, as can be seen in that figure and as noted by 
FV05, there is a metallicity offset between stars with 
and without pl anets at all masses b etween 0.7 M0 and 
1.4 M0(see also lSantos et al.l (|2004[ )V Thus, despite the 
artificial mass-metallicity correlation in our sample of 
FGK dwarfs, there still exists a clear PMC. At a given 
mass, stars with planets have higher metallicities than 
the stars without planets. 

The PMC is also apparent in the M dwarf sample. 
The low-mass stars with planets are extremely metal- 
rich compared to the full stellar sampl^H- Also apparent 
from the M dwarf sample is that there are far fewer planet 
detections, both in an absolute and fractional sense, com- 
pared to the higher-mass stellar samples. 

The far right-hand panel of Figure [2] shows that the 
metallicity offset between stars with and without planets 
is also present among the more massive subgiants, albeit 
at lower statistical significance. Like the FGK dwarfs, 

^ The metallicities of the full sample of M dwarfs were estimated 
using the photometric calibration of Johnson & Apps (2009). This 
required an extrapolation of their relationship for the stars below 
the main sequence. However, since the relationship between AMk 
and [Fe/H] is expected to be monotonic, our extrapolation will not 
affect our conclusions in this case. 



the subgiants have a artificial mass-metaUicity correla- 
tion, owing to the red cutoff oi {B — V < 1.1) used in the 
selection of the subgiants from the Hipparcos catalog. 
A much higher fraction of massive stars have detected 
planets than do the M or FGK dwarfs. 

The metallicity offsets among the stars with and with- 
out planets in the three mass-bins in Figure [5] are sug- 
gestive of a PMC that spans an order of magnitude in 
mass, from 0.2 M© to 2.0 M0. Also seen among the three 
mass bins is a steadily increasing planet occurrence rate: 
while only 3.3% of the M dwarfs have a planet, 20% of 
the retired A stars harbor one or more giant planets. 
This is strong evidence that planet occurrence correlates 
with stellar mass, separately from the effects of stellar 
metallicity. In the following sections we examine these 
trends in further detail. 

4. QUANTIFYING PLANET OCCURRENCE 
4.1. Parametric Description 

We derive a parametric relationship between stel- 
lar properties and fraction of stars with planets using 
Bayesian inference. The resulting function, while ad hoc, 
can be used to predict yields of future planet surveys, in- 
terpret the results of ongoing planet search efforts, and 
compared directly to the output of theoretical models of 
planet formation. 

Our choice of functional f orm follows from the m etal- 
Hcity analyses of FV05 and lUdrv fc Santos! ()2007| ). who 
describe the fraction of stars with planets, /, as a func- 
tion of metallicity in the form f{F) ex 10'^^, where 
F = [Fc/H]. Our parametric model also needs to ac- 
count for stellar mass. Previous observational studies 
suggest that planet oc currence should rise monotoni- 
cally with stellar mass (ILaws et al.ll2003t Ijohnson et al.l 
I2007at iLovis fc Mavori 120071 ). For the mass relationship 
we adopt a power law /(M) oc M", where M = M*/Mq. 

Since we assume that mass and metallicity produce 
separate effects, the fraction of stars with planets as a 
function of mass and metallicity can be described by 



f{M,F) = CAf"10'^^ 



(1) 



We note that there exist many possible functional 
forms for f{M) and f{F). Indeed, any monotonic func- 
tion sho uld provide an adequat e fit to our data set. For 
example. iRobinson et al.l (|2006[ ) use a logistic function to 
describe planet fraction as a function of stellar a-element 
abundance, and they note that a power law is simply an 
approximation to the low-yield tail of such a function. 
However, we have decided to use power law descriptionqf] 
due to the simplicity of the functional form and for ease 
of comparison with previous studies. 

4.2. Fitting Procedure 

For conciseness, we denote the parameters in Equa- 
tion [1] by X. The parameters can be inferred from the 
measured number of planet hosts H drawn from a larger 
sample of T targets using Bayes' theorem: 



P{X\d) o:P{d\X)P{X) 



(2) 



* Since [Fe/H] oc logNp^ , the exponential term in Equation [T] 
is a power law relationship of the number of iron atoms: f{F) oc 



where P{X \d) denotes the probability of X conditioned 
on the data d. In our analysis, the data represent a bi- 
nary result: a star does or does not have a detectable 
planet. The terms on the right of the proportionality are 
the probability of the data conditioned on the distribu- 
tion of possible X, multiplied by the prior knowledge and 
assumptions we have for the parameters. 

Each of the T target stars represents a Bernoulli trial, 
so the probability of finding a planet at a given mass 
and metallicity is given the binomial distribution. The 
probability of a detection around star i (of H total de- 
tections) is given by f{Mi,Fi). The probability of the 
jth nondetection is 1 — f{Mj,Fj). Thus, 



H 

P{X\d)^P{X)l[f{M,,F,) 

i 
T-H 

>c\{[l-}{M,,F,)] (3) 



For each detection or nondetection, our measure- 
ments of the stellar properties of each system. Mi and 
Fi, are themselves probability distributions given by 
Pohs{Mi,Fi). We approximate these pdfs as the product 
of Gaussian^ with means {Mi^Fi} and standard devia- 
tions {(JM.i,<^F,i\- The predicted planet fraction for the 
ith star can then be expressed as 



/(M„ F^)^ j j Pobs(M„ Fi)f{M, F)dMdF (4) 

For ease of calculation the products in Equation [3] 
can be rewritten as the sum of log-probabilities, or the 
marginal log-likelihood 



H 



£ = log P(d I X) a ^ log /(Af„ F,; 



T-H 



-Y,\og[l-!{M,,F,)] 

j 
\ogP{X) 



(5) 



The parameters X — {C, a, /3} are then optimized by 
maximizing C conditioned on the data. 

We perform our maximum-likelihood analysis by nu- 
merically evaluating £ on a 3-dimensional grid over in- 
tervals bounded by uniform priors on the parameters 
{C, a,/?}. In our case, the priors simply define the inte- 
gration limits on the marginal probability density func- 
tions (pdf) of the parameters, e.g. 



P{a I d) 



C,„a, 



P{X I d)dl3dC 



(6) 



c„ 



10/3^ oc N^^ 



^ Because stellar metallicities are used to select the appropriate 
stellar model grids ( "isochrones" ) for the estimate of the stellar 
mass, these two measurements are actually covariant. However, 
we find that our result is not affected by assuming independent 
Gaussians. 



TABLE 1 

Model Parameters 



Parameter 
Name 



Uniform 
Prior^ 



Median 
Value 



68.2% Confidence 
Interval 



(0.0, 3.0) 

(0.0, 3.0) 

(0.01, 0.15) 



1.0 
1.2 
0.07 



(0.70, 1.30) 

(1.0, 1.4) 
(0.060, 0.08) 



^ We used uniform priors on our parameters between the 
two limits listed in this column. 

The ranges of the uniform priors used in the analysis are 
hsted in the second column of Table [1] 

It might at first seem more appropriate to use a prior 
for /3 from the analysis of FV05, e.g. a Gaussian centered 
on /3 = 2, rather than a uniform function. However, we 
decided against this choice of prior because no confidence 
interval for /3 was reported by FV05, and we could not be 
certain that their value was truly representative of our 
data due to fundamental differences in our methodology, 
as we discuss in § [5] Similarly, no function al form for 
f{M) was reported bv lJohnson et al.1 ()2007aD . 

However, our choice of a uniform prior is not entirely 
uninformed. Based on previous studies, we felt it was safe 
to consider only monotonically increasing functions {a > 
0, /3 > 0), and for f] we chose a range that encompasses 
the value measured by FV05. 

5. RESULTS 

The best-fitting parameters and their 68.2% ("l-cr") 
confidence intervals are listed in the third and fourth 
columns of Table [TJ We estimated the confidence inter- 
vals by measuring the 15.9 and 84.1 percentile levels in 
the cumulative distributions (CDF) calculated from the 
marginal pdf of each parameter (e.g. Equation [6]) . 

The marginal joint parameter pdfs are shown in Fig- 
ure[31 The comparisons between the best-fitting relation- 
ship (Equation[T|) and the data are shown in Figure|3]and 
[5l In both figures, the histograms show the "bulk" planet 
frequency, with bin widths of 0.15 Mq and 0.1 dex, re- 
spectively. The filled circles denote the median planet 
fraction predicted by Equation [T] based on the masses 
and metallicities of the stars in each bin. The diamonds 
show the best-fitting metallicity and mass relationships, 
given by f{M,F = 0) and J{M = 1, F). 

Our Bayesian inference analysis provides two addi- 
tional assurances that stellar mass and metallicity corre- 
late separately with planet fraction. The first is the lack 
of covariance between a and /3 in Figure [31 This also 
demonstrates that our stellar sample adequately spans 
the mass-metallicity plane despite the artificial correla- 
tion between stellar parameters in part of our sample. 
The second check on our initial assumptions is seen in 
Figure ID While some of the increase in planet fraction 
as a function of stellar mass is due to a rise in aver- 
age stellar metallicity in our sample of high-mass stars 
(circles), there still exists a nearly linear increase owing 
to stellar mass alone (diamonds). Thus, there is an ap- 
proximately order-of-magnitude increase in planet occur- 
rence over the mass range spanning M dwarfs to A-type 
stars. Similarly, some of the metallicity relationship is 
due to the higher stellar masses among the metal-rich 
stars. However, there still exists a strong metallicity cor- 
relation spanning more than an order of magnitude in 




0.0 0.5 



Fig. 3. — Marginal posterior pdfs for the model parameters con- 
ditioned on the data. 

iron abundance. 

In the following section we compare our results to those 
of related studies. 

6. COMPARISONS WITH PREVIOUS WORK 
6.1. Previous metallicity studies 

FV05 studied planet occurrence as a function of metal- 
licity among the Sun-like portion of our stellar sample 
and found /3 = 2. By restricting our analysis to the 
SPOCS subset of our sample that overlaps with FV05, 
and by fitting a function of metallicity alone we find 
/3 = 1.7±0.3, which agrees with the FV05 value to within 
our 68.2% confidence interval. The significance of the dif- 
ference is reduced further if we assume the uncertainty 
in their measurement is comparable to ours. 

By fitting for both mass and metallicity we find /3 — 
1.4 ± 0.3 and a — 0.7 ± 0.4, which agree with the values 
in Table 1 measured for the full stellar sample. This pro- 
vides assurance that our analysis is not overly sensitive 
to our high-mass stellar sample, among which our detec- 
tion sensitivity is lower due to the shorter time baseline 
and fewer Doppler measurements per star. 

It is likely that this smaller (3 from our analysis of the 
Sun-like stars compared to that of FV05 is in part due to 
the different methods of fitting the planet-fraction rela- 
tionship, i.e. their least-squares fit to histogram bins ver- 
sus our Bayesian approaclCB- The other key difference 
is that we simultaneously fit to both mass and metal- 
licity. Since mass and metallicity are correlated in our 
samples some of the metallicity relationship observed by 
FV05 was due to stellar mass. This effect can also be 
seen in Figure |4l The joint mass-metallicity relation- 
ship sits above the metallicity power-law at high values 

^^ For example, FV05 performed a x^ minimization, which as- 
sumes symmetric (v(^)) error bars on their histogram bins. How- 
ever, the errors should have been binomial and asymmetric, which 
would have admitted smaller values of /3. 



of [Fe/H] since the metal-rich stars in our sample tend to 
be slightly more massive on average than the metal-poor 

st ars. 

lUdrv fc Santoi (J2007D analyzed the FV05 sample, to- 
gether with a sample of stars drawn from the CORALIE 
survey, and found /3 = 2.04 for [Fe/H] > 0. For lower 
metallicities they suggest a flat occurrence rate provides 
a better fit than the continuation of the exponential re- 
lationship to sub-Solar metallicities. We compared the 
two functional forms (exponential versus exponential- 
plus-constant) using the method of Bayesian model com- 
parison. By integrating the right-hand side of Equation[2] 
over all parameters X, one obtains the evidence, or total 
probability of the model conditioned on the data: 



P{d)^ 



P{a,/3,C\d)dadpdC 



(7) 



The ratio of evidences provides a means of quantifying 
preference in one model over another. If a model has 
evidence more than a factor of 10 greater than the alter - 
native, it is "strongly preferred" ()Kass fc Raftervl[i995[ ). 
When fitting planet-fraction as a function of metallicity 
alone, the evidence for the exponential-only model is a 
factor of 1800 higher than the exponential-plus-constant. 
Thus, our data strongly prefer a model in which the 
fraction of stars with planets continues to decrease for 
[Fe/H J < 0. 

We can take the Bayesian evidence analysis a step fur- 
ther and compare our joint fit to the planet fraction as 
a function of mass and metallicity to previous fits to 
metallicity or mass alone. We find that the evidence 
for the joint fit is a factor of 2400 larger than that of 
the metallicity-only fit, and a factor of 10^ higher than a 
mass-only fit. The planet fraction among our sample is 
therefore best described as a function of metallicity and 
stellar mass. 

6.2. Previous mass studies 

iJohnson et al.l ()2007al ) used roughly the same sample 
presented herein to measure the occurrence rate of plan- 
ets in three coarse mass bins with widths of 0.6 Mq cen- 
tered on Ml, = {0.4, 1.0, 1.6} Mq. In these three intervals 
they measured occurrence rates of 1.8 ± 1.0%, 4.2 ±0.7% 
and 8.9 ± 2.9%. After correcting for the average stel- 
lar metallicity in each bin, the fractions change slightly 
to 2.5 ± 1.2%, 3.5 ± 0.7, and a lower limit of 6.3% for 
the high-mass bin. Integrating our relationship over the 
same mass intervals yields 2.5 ± 0.9%, 6.5 ± 0.7%, and 
11 ± 2%. 

The agreement for the low-mass bin is not too surpris- 
ing since we are using the same sample of M dwarfs as 
used by Johnson et al. The disagreement for the FGK 
dwarfs is in part due to the different selection criteria for 
planet detections: we use a velocity amplitude cutoff of 
i^ > 20 m s"""^, compared to the Mp sinz > 0.8 Mjup 
used by Johnson et al. Because of this, our sample of 
planet detections includes a larger number of low-mass 
planets at short orbital periods, particularly for the Sun- 
like stellar sample. At higher stellar masses, our mea- 
sured planet fraction represents a significant refinement 
over the result of Johnson et al., which stems primarily 
from our larger sample size and higher Doppler precision 
with Keck/HIRES compared to Lick/Hamilton. 



Our revised planet fraction for the high-mass stars ap- 
pears much smaller than the recent results presented 
bv iBowler et all ()2010l ). who measured / = 26J:g% for 
1.5 < Mj^/Mq < 1.9, based on the Lick subgiants sam- 
ple. However, Bowler et al. reported the bulk occurrence 
rate, and did not attempt to correct or fit for metallicity. 
Our analysis shows that metallicity plays an important 
role in shaping the bulk occurrence rate among our sub- 
giants, which are metal-rich by -1-0.14 dex compared to 
the less massive stars. This can be seen in the highest 
mass bin in Figure HI in which the measured planet frac- 
tion is consistent with the value measured by Bowler et 
al. Simila rly, in their analysis o f the planet fraction for 
M dwarfs (jJohnson et al.|[2010bf ) noted the higher occur- 
rence for metal-rich M dwarfs, but only reported a bulk 
occurrence rate for the sample. 

6.3. Is there a planet-metallicity correlation among our 
evolved stars? 

In their analysis of the metallic i ty dis tribution of K gi- 
ants with planets. fPasquini et al.l ()2007[ ) concluded there 
was no evidence of a PMC among their evolved stars. 
iTakeda et al.l ([2008) found a similar result based on their 
sampl e of massive K giants, while .Hekker fc Melendezl 
(|2007f ) did find evidence supporting a PMC among their 
stellar sample. This somewhat contentious point has im- 
portant implications for the interpretation of the PMC 
seen among Sun-like stars. Pasquini et al. argued that 
this lack of a metallicity correlation was evidence for the 
"pollution" scenario, in which only the outer layers of 
stars with planets were metal-enriched by the infall of 
gas-dep leted planetesimals during the planet-form ation 
epoch (jGonzaled 119971 : iMurrav fc Chabove J l2002f ) . In 
this scenario, as stars evolve off of the main sequence 
their convective envelopes deepen and their polluted 
outer layers are dil uted, erasing an y "skin-deep" metal- 
licity enhancement (|Laughlinll2000D . 

In our analysis, we implicitly assume that the PMC 
holds among the evolved stars in our sample. We can 
test this assumption by restricting our analysis to M^, > 
1.4 Mq and comparing the Bayesian evidence between 
two models: planet fraction as a function of stellar mass 
alone, f{M) — CM"' (corresponding to a fiat metallicity 
distribution) versus planet fraction as a function of stellar 
mass and metallicity, f{M, F), given by Equation [TJ We 
fitted both models to the subsample of massive subgiants, 
which have deep convectiv e envelopes accordi ng to the 
Padova stellar model grids (jGirardi et al.ll200^ . We find 
that the planet fraction among these evolved stars is best 
described as a function of mass and metallicity, with an 
evidence ratio of order lO^'^. 

The extreme magnitude of the Bayes factor is driven 
primarily by the 7 subgiants with M^, > 1.4 Mq and 
[Fe/H] > +0.35 (Figure O. It is highly improbable that 
a flat metallicity distribution would result in 5 out of 7 of 
these metal- rich subgiants harboring a planet. The best 
fitting parameters are a = 1.5 ± 0.4 and (3 — 0.73 ± 0.35, 
which are lower than, yet consistent with the values we 
measure for the full stellar sample. However, the size 
and metallicity range of our sample of subgiants only 
allows us to rule out a flat metallicity relationship with 
98% confidence. At present we can say that our data are 
consistent with a PMC among our massive subgiants. 
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Fig. 4. — Planet fraction (/ = Npianats [N stars) as a function of mass for our stellar sample (gray histogram). The red filled circles show 
the planet fraction predicted by Equation[T]for the masses and metallicities of the stars in each histogram bin. Note that we use a histogram 
only for visualization purposes; the data were fitted directly without binning. The open diamonds show the best-fitting relationship between 
planet fraction and stellar mass for [Fe/H] = 0. The dashed line shows the stellar mass relationship predicted bv lKennedv &: KenvonI 12009 ) 
for Solar metallicity. 



We are not certain about the source of disagreement 
between our result and those of Pasquini et al. and 
Takeda et al. One possibility is the difference in our 
statistical methodologies. Both of those previous stud- 
ies compared the histograms of stars with and without 
planets, and as a result did not quantify their confidence 
in a PMC or lack thereof. It will be informative to apply 
the techniques outlined in § |4] to the stellar and planet 
samples of the various K-giant surveys in order to make 
a meaningful comparison with our results. 

7. SUMMARY AND DISCUSSION 

We have used a large sample of planet-search target 
stars and planet detections from the CPS to study the 
correlation between stellar properties and the occurrence 
of giant planets (K > 20 m s~i) with a <2.5 AU. We 
have derived an empirical relationship describing giant 
planet occurrence as a function of stellar mass and metal- 
licity, given by 



/(Af^, [Fe/H]) =0.07 ± 0.01 x {MJMq)^-^^^-'^ 

^ -^q1.2±0.2[Fo/H]^ 



(8) 



Our understanding of planet formation is presently 
dominated by two theories: core accretion (e.g. 
FoUack ct al..,1996) and disk instability (,Boss,1997, ). The 



core accretion model is a bottom-up process, by which 
protoplanetary cores are built up by the collisions of 
smaller planetesimals. Once the core reaches a critical 
mass of roughly 10 M^, it rapidly accretes gas from the 
surrounding disk material. The disk instability mecha- 
nism is a top-down process whereby giant planets form 
from the gravitational collapse of an unstable portion of 
the protoplanetary disk. 

Both models depend on the existence of a massive gas 
disk, a portion of which forms the bulk of the final mass 
of the Jupiter-like planet. However, since the inner gas 
disks of protoplanetary disks disperse on timescales of 3- 
5 Myr, th e process of pla.net fo rmation is a race against 
time fe.g. iPickett fc Liinll2004[ ). In this race the disk in- 
stability model holds a major advantage over the core ac- 
cretion model because, under the right conditions, plan- 
ets can form from disk collapse in a mere thousands of 
years, compared to of order Myr timescales required for 
core accretion. 

The disk instability model predicts that there should 
be no dependence on planet formati on an d physical stel- 
lar properties. The simulations of iBossI (2006) showed 
that giant planets should readily form in even low-mass 
protoplanetary disks, and that in general giant planets 
should form efficiently via disk instability independent 
of stellar mass. The disk instability model also predicts 
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Fig. 5. — Planet fraction (/ = Npianets / N stars) SlS a function of nictallicity for our stellar sample (gray histogram). The red filled circles 
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between planet fraction and stellar metallicity for M* = 1 Mq. None of the 52 stars with [Fe/H] < —0.5 harbor a giant planet. 



that planet formation should also be independent o f disk 
metalhcity (Ep ss 2002). Indeed, iCai et al.l (|2006[ ) and 
IMeru fc Bate! (i2010) show that the efficiency of disk in- 
stability to form giant planets decreases with increasing 
metallicity. 

In contrast to these predictions of the disk instability 
model, we observe a strong dependence between planet 
occurrence and the physical properties of the star. As- 
suming that the present-day mass and metallicity of a 
star reflects the conditions in its protoplanetary disk, 
then our results suggest that disk instability is not the 
primary formation mechanism for the giant planets de- 
tected by Doppler surveys. Indeed, it has long been 
recognized in that there are theoretical complications in 
forming close-in planets via disk instability ?, e.g.. How- 
ever, the m echanism may be responsible for planets in 
wide orbits (jDodson- Robinson et al.l [20091 ) . however see 
? for complications with this scenario. 

An alternative explanation for the observed planet- 
metallicity correlation is the so-called "pollution" model. 
In this scenario, planet formation actually occurs around 
stars of all metallicities, and the accretion of gas-depleted 
protoplanets onto the thin convective layers of stars 
gives rise to an enhanced stellar metallicity that i s 
actually only " skin deep" (jMurrav fc Chaboveii [20021 ). 
iPasQuini et al.l (|2007[ ) interpreted the flat metallicity dis- 
tribution of K giants with planets as evidence for such 



an effect, since the deepening convective envelopes of 
evolved stars should dilute any metallicity enhancement 
of its outer layers. However, we find that the planet frac- 
tion among our massive subgiants is described well by a 
model with a monotonic rise as a function of both mass 
and metallicity (Section [5]) . 

While much attention has been given to evolved, 
intermediate-mass star s in the in vestigation of the pol- 
lution paradigm (Laug hlin fc Ada ms 1997, FV05), stars 
at the other end of the stellar mass scale provide an- 
other proving ground. M dwarfs have deep convective en- 
velopes over their entire lifetimes, and stars with masses 
below 0.4 M© are expected to be completely convective, 
at least in the absence of s trong magnetic activity (e.g. 
iMuUan fc MacDonaldl [2007). 

The evidence for a PMC among the subgiants, to- 
gether with a strong PMC seen among the M dwarfs, 
are highly suggestive that the present-day metallicities 
of stars are representative of the compositions of their 
disks during the planet-formation era. One compelling 
explanation for both the observed PMC and stellar-mass 
correlation is that the surface density of solids is a key 
factor in the planet f ormation process ([Laughlin et al.l 
[200l [IdalFLg [2005 [Robinson et al.ll2006l) . If so. both 
higher stellar (disk) metallicity and higher stellar (disk) 
mass can generate the requisite surface density for planet 
formation. 
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The relationship between planet formation efficiency 
and stellar mass/metallicity has been previ ously studied 
in the context the core accretion paradigm flkoma ct alj 
2000; Korn et et_alT 2005. 2006; Kennedy & Kcnvon 200l 
Thommcs et al. 2008; Krctkc ct al. 2009; Johansc n et aLl 
2009; Mordasini ct al. 2009; Dodson-Robinson et alJ 
20091 ). Kennedy fc KenvonI (|2008l ) modeled the evolu- 
tion of the temperature profile at the disk midplanes of 
stars of various masses, and studied the width of the ra- 
dial region of disks in which protoplanetary cores form 
most efficiently. They found that the disks around A- 
type stars on their descent to the main sequence have 
very broad formation regions, and their models pre- 
dicted a positive correlation between stellar mass and 
giant planet formation efficiency. Their prediction for 
the planet fraction is shown as a dashed line in Figure |4l 
which we have approximated using the polynomial rela- 
tionship /(M) = -0.03633 -H0.0138A/^-0.0060Af 2. The 
agreement between theory and observation is striking. 

The interplay between the mass and metallicity of pro- 
toplanetary di sks is also apparent in the core accretion 
simulations of iThommes et aTl (J2008D . In their analysis, 
they simulated disks with a wide variety of masses, vis- 
cosities and metallicities. Their models produce gas gi- 
ants most effectively in disks with a combination of high 
masses (Mdisk ^ 0.04 Mq) and low viscosities. In their 
simulations of the effects of metallicity, they found that 
gas giants can form in Solar-composition disks only if the 
disk masses exceed ~ 0.06 Mq, or twice the minimum- 
mass Solar nebula. This mass threshold decreases to 
« 0.03 Mq for disks with [Fe/H] = 0.25. Thus, Thommes 
et al. showed that there can be a trade-off between the 
mass and metallicity of a protoplanetary disk in forming 
giant planets. In the core accretion paradigm, M dwarfs 
can form giant planets, but only if they have high metal- 
licities. Similarly, even low-metallicity A stars can form 
mas sive planets, owing to their m ore massive disks (see 
also iDodson-Robinson et al.l I2009D . 

Theories of planet formation have progressed in large 
leaps thanks largely to the rapidly growing sample 
of exoplanets discovered around other stars. Suc- 
cessful, generalized theories of the origins of plane- 
tary systems must account for the observed correla- 
tions between planet occurrence and stellar properties. 
The findings presented herein are the result of more 
than 15 years of high-precision Doppler monitoring of 
nearby stars. Additional information will soon pour in 
from other surveys using techniques such as microlens- 



ing (e .g. iDong et al.l 120091; iGould et al]|2010D. a. strom- 
etrv (iBoss et all l2009D~tra nsits (B orucki et all 120041; 
iBarge et al."2008'; 'Irwin et al.ii20 09) and direct imaging 
jClaudi ct al. 2006; Arti gau et al . 2008; M acintosh et al.l 
[2008). 

Our results have important implications for the target 
selection of these future planet surveys. We find that 
A-type stars harbor planets at an elevated rate com- 
pared to less massive stars. However, it is not obvious 
whether A type stars make the most promising targets 
for other types of surveys. Where massive stars per- 
haps hold the most promise is for direct imaging surveys. 
The first two planetary systems imag ed around normal 
stars were both yo ung A-type dwarfs (jKalas et al.l 120081 ; 
iMarois et al.ll2008l ). Are A dwarfs the ideal direct imag- 
ing targets? The answer to this question will rely on the 
results presented herein, along with a careful consider- 
ation of the mass, metallicity, luminosity and age dis- 
tributions of nearby stars; and the orbital and physical 
properties of planets as a function of stellar mass. This 
issue will be addressed in a companion paper (Crepp & 
Johnson 2010). 
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